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Abstract. We investigate the evolution of spherically symmetric supernova remnants in which mass loading takes 
place due to conductively driven evaporation of embedded clouds. Numerical simulations reveal significant dif- 
ferences between the evolution of conductively mass loaded and the ablatively mass loaded remnants studied in 
Paper I. A main difference is the way in which conductive mass loading is extinguished at fairly early times, once 
the interior temperature of the remnant falls below ~ 10 7 K. Thus, at late times remnants that ablatively mass 
load are dominated by loaded mass and thermal energy, while those that conductively mass load are dominated 
by swept-up mass and kinetic energy. Simple approximations to the remnant evolution, complementary to those 
in Paper I, are given. 
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1. Introduction 

The evolution and properties of supernova remnants are 
fundamental to many areas of astrophysical research. Self- 
similar solutions for various stages of the evolution have 
been obtained for various assumptions, as have unified so- 
lutions which link these together (see Cioffi, McKee & 
Bertschinger 113551 Truelove & McKee IT555)l . The self- 
similar and unified solutions are complemented by a wide 
range of numerical investigations. Amongst the problems 
examined numerically are explosions in plane-stratified 
media (e.g. Falle & Garlick Arthur & Falle ITWT1 

1993), the evolution of SNRs inside pre-existing wind- 
driven cavities (Franco et al.[1991), the influence of hy- 
drodynamic instabilities (e.g. Chevalier & Blondin 1995), 
and applications to the broad emission line regions in ac- 
tive galactic nuclei (e.g. Terlevich et al. 119921 Pittard et 
al. EOcfl)) . 

In most investigations of supernova remnant evolu- 
tion, the surrounding medium has been assumed to have 
a smooth density distribution. However, the interstel- 
lar medium is known to be multi-phase (e.g. McKee & 
Ostriker 1977), and there is clear observational evidence of 
engulfed clouds within the supernova remnant N63A (Chu 
et al. 119991 Warren et al. l2002(l . As it is widely known, the 
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injection of mass from such cool clouds into the interior of 
a supernova remnant will affect its behaviour and struc- 
ture. 

If the clouds are assumed to be continuously dis- 
tributed, similarity solutions describing the evolution of 
adiabatic mass loaded supernova remnants can be derived 
for a specific functional dependence of the mass injec- 
tion. Previous work has focussed on cases in which one 
mass injection process is dominant over others. McKee & 
Ostriker {1377) 1. Chieze & Lazareff (PHD, and White & 
Long l|1991[) each investigated the evolution under the as- 
sumption that cloud destruction is due to conductive evap- 
oration. Alternatively, Dyson & Hartquist (|1987H treated 
hydrodynamic ablation as the dominant process. 

A small number of papers based on numerical simula- 
tions of mass loaded supernova remnants also exist in the 
literature. Cowie, McKee & Ostriker l|1981(l included the 
dynamics of the clouds and found that warm clouds are 
swept towards the shock front and are rapidly destroyed, 
while cold clouds are more evenly distributed and have 
longer lifetimes. When energy losses become important, 
this behaviour leads to the highest densities in the rem- 
nant occuring over the outer half radius. In contrast, when 
there is no mass loading, a thin dense shell forms at the 
forward shock. Arthur & Henney ©96 ) studied the effects 
of mass loading by hydrodynamic ablation on supernova 
remnants evolving inside cavities evacuated by the stellar 
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winds of the progenitor stars. They showed that the extra 
mass injected by embedded clumps was capable of produc- 
ing the excess soft X-ray emission seen in some bubbles in 
the Large Magellanic Cloud. 

The range of a supernova remnant is relevant to a num- 
ber of important areas of astronomy. It can affect the effi- 
ciency of sequential star formation and the global dynam- 
ics and structures of starburst superwinds where many 
remnants overlap. The superwind in the starburst galaxy 
M82 must be mass loaded to account for the observed X- 
ray emission (Suchkov et al-HHU- Since the range and 
radiative energy losses of a remnant are affected by mass 
loading, it is therefore desirable to have approximations 
which describe remnant evolution and range in clumpy 
media. 

Analogous approximations already exist for smooth 
media (cf. Truelove & McKee 1999]and references therein). 
In a recent paper (Dyson, Arthur & Hartquist 2002 , here- 
after Paper I), the first steps towards this goal were taken 
with the derivation of range approximations for cases in 
which mass injection occurs at a constant rate, q, or at a 
rate which depends on the flow Mach number relative to 
the clump (Hartquist et al. H986[) . A mass loading rate of 
gM 4 / 3 was used for subsonic flow, where M is the Mach 
number of the flow. A constant mass loading rate q was 
used where the flow is supersonic. This prescription sim- 
ulates mass stripping from the clump by hydrodynamic 
ablation. 

In this companion paper we consider mass injection 
from clouds into the hot remnant interior by conductive 
evaporation (cf. Cowie & McKee 119771 McKee & Cowie 
EHEJ. This process is potentially important in any astro- 
physical system where a hot tenuous phase coexists with 
a colder denser phase (see McKee & Cowie I1U77I for some 
examples), though it can be suppressed by magnetic fields. 
In many systems (including supernova remnants) the hot 
phase flows past the clouds, and they will be subject also 
to ablation. Unfortunately, we currently lack even basic 
models of this combined interaction, so it is difficult to 
know whether one process dominates the other, or if one 
process limits the other. For example, one could imagine 
a scenario whereby ablation drives the initial mass loss, 
increasing the effective surface area of the cloud (e.g. by 
forming a tail), at which point conduction takes over and 
leads to efficient mixing. In evolving objects such as super- 
nova remnants, one process may dominate at early evolu- 
tionary stages, while the other dominates at later stages. 
Given these potential complications and our lack of knowl- 
edge about them, as a first approach we treat each process 
individually. 

The outline of our paper is as follows. In Sec.[21we out- 
line our assumptions and our numerical code; in Sec.|21we 
present some of our results; in Sec.0]we present simple ap- 
proximations to these results; and in Sec. we summarize 
and conclude. 



2. The Calculations 

As initial conditions, we take freely-expanding cold ejecta 
at constant density (typically 10 4 times the ambient den- 
sity, no), and with a linear velocity profile. The maxi- 
mum ejecta speed is 4000 kms" 1 , which together with 
the radius of the ejecta defines the initial time. We per- 
formed extensive tests with different density ratios to con- 
firm that our chosen ratio is large enough not to influence 
the resulting evolution. The calculations were performed 
using an adaptive grid hydrodynamic code (see e.g. Falle 
& Komissarov 119961 I1998[l . and we have further verified 
that the resolution employed is high enough to give robust 
results for the global quantities considered in this paper. 

All calculations were performed in spherical symme- 
try, with an ejecta mass of 10 Mq and a kinetic energy 
of 10 51 ergs (the effect of different progenitor masses and 
explosion energies is commented on in Sec. EJ. The as- 
sumption that the energy is entirely kinetic for the initial 
conditions is somewhat different to that in Paper I, but it 
does not affect the resulting evolution (see also Gull [1973 
and Cowie et al.[l981). We modelled the intercloud am- 
bient medium as a warm plasma of temperature 10 4 K. 
Cooling appropriate for an optically thin plasma of so- 
lar composition and in collisional ionization equilibrium 
was included. We assumed that cooling below 10 4 K is 
balanced by photoionization heating from the diffuse field 
and for that reason impose a temperature floor of 10 4 K. 
For simplicity we did not include magnetic fields or shock 
precursors, and also did not consider conduction between 
the hot remnant interior and cold neighbouring gas. 

The rate of mass evaporation from embedded clouds is 
dependent on many factors (cf. Cowie & McKee H9771) . in- 
cluding the temperature of the hot phase, the cloud radius, 
and the presence of magnetic fields. However, motions of 
the clouds are generally not large enough to change the 
evaporation rates, and Cowie & Songalia (1977) noted 
that nonspherical clouds may be treated in an approxi- 
mate way by adopting half the largest dimension as the 
radius of the cloud. The evaporative mass-loss rate from 
a single cloud (Cowie & McKee IBF7?j) is 

m = 47ra 2 nT 1/2 (fc^m H ) 1/2 0F(cr o ), (1) 

where a is the cloud radius, n and T are the density and 
temperature of the surrounding hot phase (i.e. the interior 
of the remnant), fimn is the mean mass per particle, and 
<fi is an efficiency factor which encapsulates the magnetic 
field strength, its configuration, the cloud geometry, etc. 
(0 = corresponds to no evaporation, (f> « 1 to evapora- 
tion in the absence of magnetic fields - e.g. 0=1.1 for a 
plasma with equal electron and ion temperatures and cos- 
mic abundances). F(<jq) is a function defined in Eqs. 61 
and 62 of Cowie & McKee (|1977[> . In the classical limit, 
F(ao) — 2<Jo, where the saturation parameter, cto, is de- 
fined as 

ao = ( 1.54 x10V (2) 



Pittard et al. : Evolution of Mass Loaded SNR II. 



3 



Here 



a /(I pc) is the radius of the cloud in parsecs. 



The onset of saturation occurs when ao is of order unity, 
varying between 1.95 (<fi -C 1) and 1.08 (0 = 1) for the 
examples given in Cowie & McKee (|1977[l . In the classical 
limit, rh cx T 5 / 2 since F(<jq) cx ao (cf. Fig. 4 in Cowie 
& McKee lP377|) . If saturated, F(a ) oc a^, where /3 = 
(1 + M 2 )/(6 + M 2 ), and M s (the Mach number of the 
saturated zone) is related to by 



M„(l + -A/ 2 ) = 20. 



(3) 



Since is unlikely to exceed unity by a substantial 
amount, and has a lower limit of zero, the likely range of 
/3 is 1/6 < (3 < 3/8. Hence in the saturated limit m cx T a 
where 5/6 < a < 5/4. 

Since the onset of saturation is dependent on the radius 
of the cloud for a specified hot phase, larger clouds will 
tend to evaporate in the classical limit, while the evapo- 
ration of mass from smaller clouds will tend towards satu- 
ration. As the temperature and density of the remnant 
interior evolves, the radius of clouds which are just at 
the onset of saturation will also change. Moreover, we ex- 
pect the distribution of clouds to evolve with time, as the 
timescale for the smaller clouds to completely evaporate is 
shorter than the equivalent timescale for larger clouds and 
in many cases the evolutionary timescale of the remnant. 
Further, we may expect some differences in the physical 
make-up of clouds of differing radii (e.g. the distribution 
of cloud mass in the cold core and its warmer skin - cf. 
McKee & Ostriker ll977fl . Finally, a full treatment also re- 
quires specifying the number of clouds per unit volume, 
which is likely to be dependent on the local environment 
(e.g. the cloud spectrum inside a starburst superwind is 
likely to be somewhat different to that in our local ISM). 
McKee & Ostriker l|1977f) argued that mass loading in the 
Galactic ISM would be dominated by the smaller clouds. 

Addressing all of these issues is beyond the scope of the 
present paper, so we choose instead the following simplifi- 
cation. We assume that we can specify a temperature for 
the onset of saturation. T sat , which is an average over both 
the cloud distribution and time. We adopt T sat = 10 7 K 
for all of our calculations. Given the uncertainties intro- 
duced by this assumption (and those in the original work 
of Cowie & McKee H977fl it seems unrealistic to specify 
a complicated dependence for the mass evaporation rate 
per unit volume above T — Tsat, and we therefore assume 
that it is independent of T. Tests with other values of T sat 
have shown that our results only become sensitive to the 
value of T sat when T sat < 10 7 y / n7j'K (no (cm" 3 ) is the in- 
tercloud ambient number density) . Hence our results with 
no < 1 cm~ 3 are robust (since a lower value of T sa t im- 
plies that most of the mass loading is from a population 
of tiny clouds (a pc < 0.42/0) whereas the smallest cloud 
radii in models of the ISM by McKee & Ostriker <|1977[) are 
~ 0.38 pc). On the other hand, for no > lcm~ 3 , our re- 
sults are dependent on our chosen value of T sat . However, 
as T sat will depend on details such as clump lifetimes and 



number distributions, which are problem specific parame- 
ters, we simply note that our results begin to suffer from 
loss of generality in this regime. 

Our final assumptions are: i) that the intercloud spac- 
ing is small compared to the scale over which the prop- 
erties of the remnant varies (so that a continuous mass 
source term can be used), ii) that injection takes place at 
zero velocity relative to the global flow, and iii) that the 
injected gas has zero internal energy. The full set of gas 
dynamic equations is then 



dp 1 d 2 

Ot r z Or 

dpu 1 d 2 2 2p 

Ot r z Or r 



oe 1 d r 2 . .-, 

di + ^d-ri" <e+p)]=-L, 



(4) 
(5) 
(6) 



where the symbols p, u, and p have their usual meanings, 
e = P/(l ~ 1) + 0.5/ou 2 is the total energy, and L is the 
radiative cooling term, q is the mass loading rate per unit 
volume and time. 

As in Paper I, we define a fiducial mass loading rate 



6po 

5i a f ' 



(7) 



which gives the mass flux through the blast wave divided 
by the remnant volume at the onset of thin shell formation 
for Sedov expansion. Here, p n is the intercloud ambient 
density and t s { is the timescale for thin shell formation for 
a remnant propagating into a uniform ambient medium 
given by Franco et al. (1994), 



(8) 



t s{ = 2.87 x 10 4 El{ 14 ng 4/7 yrs, 



where 10 51 i?5i is the explosion energy. In these units, qo = 
1.32 x lO- 36 ^ 3714 ^ 177 (with mean mass per particle, 
p = 0.6). 

Following the above discussion, we define the mass 
loading rate due to conduction as 



q = 316/'g 
q= fq (T/T 6 )^ 
q= 



(T > 10 7 K), 

(10 5 K < T < 10 7 K), 

(T < 10 s K), 



(9) 



where T 6 = 10 6 K (e.g. Cowie et al. lP3ST)l . Thin shell for- 
mation occurs when the temperature immediately behind 
the leading shock is about 10 6 K. At low temperatures, 
mass loading is so weak that the precise temperature cut- 
off is unimportant. Defining /' as we have in Eq.[!5] ensures 
that mass loading starts to influence the dynamics and 
evolution of the remnant if /' ~ 1. To be more specific, 
low values of /' mean that mass loading will not dominate 
in the remnant before the time of onset of thin shell for- 
mation, while /' > 1 indicates that the evaporated mass 
becomes important in the remnant before t s { is reached. 
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Readers who refer to Paper I will note that in that paper 
we defined a parameter / which is somewhat analogous to 
/'. In an ablatively mass loaded remnant the ablated mass 
is comparable to the swept- up mass at time i s f if / ~ 1. 
However, under the assumptions we have made, in very 
young remnants in which the temperature is everywhere 
above 10 7 K, / ~ 316/' for the mass loading rate in an 
ablatively loaded remnant to be comparable to that in a 
conductively loaded remnant for which Eq and no are the 
same. 

/' is physically related to the number density of clouds 
and the rate at which mass is evaporated from these. 
However, from Eq. 13 of Cowie et al. (198TJ, and since 

ci-3/14 11/7 n i , -i 

<7o oc E 51 n , we find that 

/' cx c/m pc N c E 3 5 i 14 n- 11/7 , (10) 

where N c is the number density of clouds. Thus we might 
expect real remnants to be better described by small /' 
when they evolve in environments with high ambient den- 
sities, and larger /' when they evolve in environments with 
low ambient densities. We note that if the distribution of 
clouds evolve with time, /' itself may evolve with time, 
though such an effect is not investigated in the models 
presented in this paper. 

We note here, that differences exist between this work 
and Paper I. In the calculations presented here, mass load- 
ing is permitted throughout the remnant, including the 
ejecta material, whereas in Paper I it is switched off in 
the ejecta. This causes differences when most of the mass 
loading occurs before the remnant has swept up much in- 
tercloud mass (i.e. when /' and no are large). A further 
difference is the assumed temperature of the ambient gas. 
In Paper I, a temperature of 100 K was used, whereas here 
we use 10 4 K. This will have negligible consequences un- 
til near the merger of the remnant with the surrounding 
gas. We note that in the specific case of superwind gen- 
eration, merger may take place with gas at even higher 
temperatures. 

We have computed models with /' = 
0,0.316,3.16,31.6,316, and ambient densities 
n = 0.01,0.1,1,10,100 cm~ 3 . Although there is 
an indication that high values of /' will be favoured 
when remnants evolve in surroundings with low ambient 
densities, /' also scales with the number density of clouds, 
and thus could be completely independent of the density 
of the intercloud medium. Hence we explored (/', no) 
parameter space. 

The remnants are evolved until the blast front degen- 
erates into an acoustic wave. We adopted the same defini- 
tions for the ends of the free expansion and Quasi-Sedov- 
Taylor phases as given in Paper I. These are respectively 
when the internal energy reaches 60% of the initial energy 
(FE) and when 50% of the initial energy has been radi- 
ated (QST). For /' < 3.16, the mass loading dominates 
in the remnant at times before the onset of thin shell for- 
mation, t s t, though is always decreasing in importance at 
this point (i.e. the ratio of loaded mass to swept-up mass 



is falling). For /' > 3.16 the mass loading also remains 
dominant at the time of the end of the QST stage, iQST- 
Finally, as in Paper I, we note the caveat that we are con- 
sidering a purely general case and that, in practice, /' is 
determined by the actual physical situation being stud- 
ied. We therefore ignore details such as clump lifetimes 
and spatial distributions, which only introduce additional, 
problem specific, parameters. 



3. Results 

We note that all of the results presented in this paper 
and Paper I are for a particular set of explosion param- 
eters, namely the canonical values E = 10 51 erg and 
M — 10 Mq. Where results for hydrodynamic ablation 
or constant mass loading are presented these have been 
calculated with the same code as used in the conductive 
mass loading cases but for the mass loading prescription 
given in Paper I. 

3.1. Remnant Expansion Speeds and Ranges 

In Fig.^ we show the remnant expansion speed and radius 
as functions of time for calculations with zero mass load- 
ing (/' = 0), substantial mass loading (/' = 3.16), and 
high mass loading (/' = 31.6), and for ambient densities 
of no = 0.01 cm~ 3 and n Q = 100 cm~ 3 . Mass loading 
causes a reduction in expansion velocity and range at in- 
termediate ages, but there is convergence in speeds and 
ranges at later times, regardless of the value of /' or the 
ambient density. This is in contrast to remnants which 
mass load through ablation (Paper I) where the reduc- 
tions in expansion speed and remnant range persist right 
through to the point of merger with the ambient medium. 

This dynamical behaviour can be understood as fol- 
lows. At early times, the highest temperatures are in the 
shocked ISM material so mass loading is highest there. By 
conservation of momentum, the addition of mass to this 
gas reduces its velocity and also has the effect of increasing 
the pressure (by conservation of energy). This causes the 
shocked material to brake harder and leads to an over- 
all reduction in the expansion velocity of the remnant. 
Once the ejecta have completely thermalized, the high- 
est temperatures are in the shocked ejecta material and 
most of the mass is added here, increasing the density 
and pressure in this gas. At this stage the rate of decrease 
of the expansion velocity is less in the conductively mass 
loaded remnants than in the zero mass loading case be- 
cause of the higher pressure in the remnant interior. As the 
blast wave decelerates, the postshock temperatures and 
hence the conductive mass loading rates decrease, so dur- 
ing this stage the swept-up mass becomes more important. 
Eventually, the swept-up mass dominates the dynamics in 
the postshock region and the remnant evolution tends to- 
wards the zero mass loading case. 
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3.2. Mass Loading History 

Fig. El shows the ratio of injected to swept-up mass as a 
function of time for various values of /' and ambient den- 
sities ri{). The top panels show the evolution when mass 
is injected through hydrodynamic ablation, while the bot- 
tom panels show the results obtained when mass injection 
occurs by conductive evaporation. It is immediately obvi- 
ous from this figure that these mass loading prescriptions 
yield vastly different behaviour. 

When mass loading occurs through ablation, the in- 
jected mass becomes increasingly dominant over the swept 
up mass as time progresses. The limit of the ablative 
mass loading is a constant mass loading rate q (corre- 
sponding to supersonic flow throughout the remnant), so 
that the injected mass increases as Mi oa d ~ 4:irR(t) 3 q/3, 
whereas the rate of increase of the swept-up mass is 
Mswept = <inR(t) 2 p dR(t)/dt, so that 

^load R{t)q nu 

—■ oc oc / t/t s f, (11) 

M swept R(t)p 

where we have used Eq. and approximated the rem- 
nant expansion speed as R(t) oc R(t)/t. Hence, as the 
ratio of the rates increases approximately linearly with 
time, we expect that the ratio of the injected to swept-up 
mass will also increase approximately linearly with time. 
This is indeed the behaviour that we see in Fig. [21 (a,b). 
Furthermore, for a given /', we expect that the ratio of 
the loaded mass (Mi oa d) to the swept -up mass (-M swe pt;) as 
a function of time (in units of t s {) will be independent of 
no, as is indeed seen (Fig. 121(b)). With reference to Eg. 1101 
however, "mathematical" independence of no (at constant 
/') does not imply "physical" independence (i.e. keeping 
constant the other physically meaningful parameters, like 
0, a pc , N c , E bl ). 

When mass loading occurs via conductively-driven 
evaporation, on the other hand, the behaviour of 
-Mioad/^swept as a function of t is quite different (cf. Fig. [21 
a,c). Although the injected mass exceeds the swept-up 
mass at intermediate times, as the remnant ages the rate 
of mass loading becomes negligible, and the swept-up mass 
gradually dominates. In Fig. [21 (c,d), the swept-up mass 
is, in all cases, greater than the injected mass at the point 
when the remnant merges with the ambient medium. 

This qualitative difference in the time-dependent be- 
haviour of the mass loading between the conductively 
driven evaporation and ablation cases accounts for the dif- 
ferences in the evolutionary behaviour displayed in Fig. 
of this paper and the evolution of the remnants discussed 
in Paper I. For the conductive case, the mass loading 
is extinguished at relatively early times and hence does 
not have much effect over the later stages of the remnant 
evolution. Once mass loading is "switched off" the ratio 
-Mioad /Mswept depends solely on M swept and thus is not 
independent of no. 



3.3. Mass and Energy Fractions 

Fig. El shows the time variation of the mass and energy 
fractions for no = 0.01 cm~ 3 and /' = 0, 3.16 and 31.6. 
Although our initial conditions specify the remnant en- 
ergy as being entirely kinetic at early times, a significant 
fraction is quickly thermalized. 

Values of /' = 3.16 and 31.6 correspond to / values of 
10 4 and 10 5 , respectively. These are much higher than the 
/ values discussed in Paper I. The reason for this choice 
is to approximately match the loaded to swept-up mass 
ratio in the two different mass loading scenarios at the 
onset of thin-shell formation, t a f, i.e., Mi oac i /Af SWGpt ~ 4 
in both cases. Since conductive evaporation is saturated 
only at early times in the remnant evolution (cf. Fig. [3J; 
and e) and becomes negligible once average temperatures 
fall below 10 6 K, while hydrodynamic ablation becomes 
more important as the remnant's volume increases, obvi- 
ously far higher values of / are necessary in the conductive 
evaporation case to give Mi oa d/M swcpt « 4 at i s f. 

With /' = (no mass loading) the thermal energy 
fraction peaks at approximately 0.72 i.e. the value for a 
Sedov- Taylor remnant expanding into a uniform medium 
(cf. Fig. [3h>) • Once the age of the remnant is comparable 
to i s f, the remnant begins to radiate its thermal energy, 
and the thermal energy fraction falls as the kinetic energy 
fraction rises. Mass injection speeds up these processes, 
and the thermal energy fraction first overshoots (see S I3.1|) , 
then undershoots the Sedov- Taylor value (Figs EJi,f ) . The 
decrease in the thermal energy fraction occurs not long 
after the end of the FE stage. This is long before radiative 
losses become significant, and the remnants are at this 
stage cooling by adiabatic expansion, such that thermal 
energy is converted into kinetic energy. This is despite 
the countering effect of mass loading, which at any given 
time tends to increase the thermal energy fraction relative 
to the kinetic energy fraction. However, we note that the 
thermal energy fraction begins to decline when the mass 
fraction of evaporated material starts to drop off and thus 
is simply a manifestation of the fact that mass loading is 
ceasing to be an important process in these remnants long 
before the onset of radiative cooling. In contrast, for the 
ablation cases discussed in Paper I, it was found that since 
mass loading becomes more important for the remnant as 
time goes on, the thermal energy fraction continues to 
increase right up until the end of the QST stage. 

At very late times the energy in a SNR with mass 
loading due to conductive evaporation is predominantly 
kinetic, whereas in Paper I it was found that in the fi- 
nal stages of a SNR with mass loading due to ablation 
the continued mass loading ensures that thermal energy 
dominates. 

3.4. Evolution of Physical Quantities 

In Fig. 21 we show the density, pressure, temperature, and 
velocity distributions at specific times just before and after 
the predicted onset of thin-shell formation for models with 
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n = 0.01 and /' = and /' = 3.16. The "structure" at 
small radii (where the ejected mass is) seen in the density 
and temperature plots is a consequence of the initial con- 
ditions and imposed spherical symmetry (which requires 
a reflection condition at r = 0) but does not significantly 
affect global properties (see e.g. CiofH et al. ll988(l . 

For /' = the distributions of density etc. in the exter- 
nal part of the remnant (where the swept up mass is) cor- 
respond to the standard Sedov- Taylor solution (although, 
as in Cioffi et al. 1988, at high ambient densities the rem- 
nant becomes radiative before fully entering this phase). 
The thin-shell formation time, i s f, for this ambient den- 
sity is 4 x 10 5 yr (from Eq. |HJ. Inspection of Fig. 0] (in 
particular, the pressure, which drops markedly in the re- 
gion where strong radiative cooling starts) shows that our 
simulations are consistent with this value. 

The addition of mass through the conductively driven 
evaporation of embedded clouds significantly alters the 
properties of the remnant, most obviously causing the in- 
terior density of the remnant to increase with respect to 
the no mass loading case. At early times, mass loading oc- 
curs mainly in the centre of the remnant (in the shocked 
ejecta region). Once temperatures here fall below 10 6 K 
(due to adiabatic expansion) , mass loading occurs mainly 
in the region of hot post-shock gas behind the blast wave, 
creating a "thick-shell" morphology (cf. Fig. Et) , resem- 
bling that of Cowie et al. (1981) and Dyson & Hartquist 
lP57|> . 

3.5. Volume- Averaged Mass Loading Properties 

In Fig. 0(a) we show the remnant-averaged mass-injection 
rate (rate of mass evaporation divided by the volume of 
the remnant) as a function of time. At early times, the en- 
tire remnant is above T = 10 7 K, hence conductive evap- 
oration is saturated and the remnant-averaged mass load- 
ing rate is approximately constant as the remnant expands 
(the small rise is caused by the increasing fractional vol- 
ume of the hot shocked gas in the remnant). Once the av- 
erage remnant temperature drops below 10 K the mass 
evaporation rate is no longer saturated and p, the rate 
of mass evaporation per unit volume, varies with time as 
p oc t~ 2 . This matches the behaviour found by Chieze & 
Lazareff (19SU, who adopt p oc T 5 / 2 , for heavily mass 
loaded remnants. This is a steeper dependence than the 
p oc t^ 1 dependence found by White & Long (1991), which 
is a consequence of their particular description for mass 
loading (pocT 5 / 6 ). 

From Fig. and Fig. we see that the majority of 
the mass evaporation from the embedded clouds occurs 
after the average temperature of the remnant drops below 
10 7 K. Once the average remnant temperature drops be- 
low 10 5 K mass loading is effectively "switched off" , the 
remnant averaged mass injection rate drops sharply, leav- 
ing the power-law dependence on time, and the quantity 
of evaporated mass remains constant in the remnant from 
this time onwards. Interestingly, this quantity does not 



Table 1. The total mass, Mfe(M q ), in the remnant at 
the end of the free expansion stage for conductive mass 
loading. Apart from models with /' = 0.316, the mass is 
relatively invariant with no, being ps 56 Mq for /' = 
and fa 37 M for /' > 3.16. E = 10 51 erg, M = 10 M Q . 









/' 






no 


0.0 


0.316 


3.16 


31.6 


316 


0.01 


56.4 


45.6 


38.1 


36.2 


36.4 


0.10 


56.4 


43.3 


37.4 


36.9 


37.6 


1.00 


56.0 


42.4 


37.6 


36.2 


36.2 


10.0 


55.6 


41.6 


36.9 


37.6 


37.6 


100.00 


56.0 


39.3 


37.6 


37.0 





depend linearly on the mass loading factor /' (Fig. |SJd). 
Roughly half the evaporated mass is loaded after the av- 
erage temperature has dropped below 10 6 K for the model 
with no = 0.01 cm~ 3 and /' = 3.16. Hence our results 
should not be strongly affected by the precise value of 
the temperature T sat at which conductive evaporation be- 
comes saturated, or by the initial conditions that we spec- 
ify- 

3.6. End of the QST Stage 

In Fig. the ratio of the time at which the total en- 
ergy of the remnant falls to one half of its initial energy, 
ii/2(/', no), to the shell formation timescale, i s f(no), de- 
fined in Eq. |H1 is shown as a function of the mass loading 
parameter /' for the five different values of the ambient 
density used. This figure shows far more variation with no 
than the equivalent figure in Paper I, where the curve for 
the different values of the mass loading parameter lie on 
top of each other. In the conductively driven mass load- 
ing case the mass loading rate is sensitive to the thermal 
structure throughout the remnant, and radiative cooling, 
when it begins, occurs in a larger volume of gas (because of 
the "thick-shell" morphology). The radiative cooling rate 
thus depends on the density and temperature of the gas 
in this extended region, which have a non-simple relation 
to the mass loading factor /'. In contrast, in the ablation 
cases discussed in Paper I, mass loading occurs in the nar- 
row region just behind the blast wave, and so the amount 
of radiative cooling in this zone (which ultimately defines 
£qst) depends directly on the factor /. 

4. Discussion 

4.1. The Free Expansion (EE) Phase 

Table gives the total mass in the remnant at the end 
of the FE phase, Mfe, for remnants with an initial en- 
ergy, E — 10 51 erg, and ejecta mass, M = 10 M Q . 
When there is no mass loading the total mass is in- 
dependent of the ambient density (being approximately 
56 M Q , cf. Paper I), though with intermediate rates of 
mass loading (e.g. /' ps 0.316) Mfe becomes dependent 
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Fig. 1. Expansion velocity (kms -1 ) and radius (pc) versus time (yrs) for remnants (E = 10 51 erg, M = 10 Mg) 
expanding into an ambient density of no — 0.01 cm -3 and no = 100 cm~ 3 . The temperature of the ambient medium is 
10 4 K, and its sound speed is 15kms _1 . The solid line has /' = (no mass loading), while the dotted and dashed lines 
have /' = 3.16 (substantial mass loading) and /' = 31.6 (high mass loading) respectively. The diamonds, triangles and 
squares indicate the end of the FE and QST stages for /' = 0, 3.16, 31.6 respectively. 



Table 2. The radius (pc) of a conductively mass loaded 
remnant at the end of the Quasi-Sedov- Taylor stage. E = 
10 51 erg, M = 10 M©. 









/' 






n 


0.0 


0.316 


3.16 


31.6 


316 


0.01 


214 


201 


167 


117 


50.7 


0.10 


85.8 


72.6 


56.4 


23.9 


7.97 


1.00 


32.8 


21.8 


11.7 


4.45 


2.00 


10.0 


12.4 


5.67 


2.77 


1.32 


0.67 


100.00 


4.35 


1.74 


0.89 


0.45 


0.23 



Mfe ~ 37 M©, and is insensitive to both /' and no. The 
fact that Mfe is dependent on /' for a given no is con- 
sistent with Figs. El (b,d,f), where it can be seen that the 
remnant energy thermalizes more rapidly with increasing 
/'. This behaviour arises from the fact that mass loading 
reduces the rate at which the remnant expands at early 
times, as it causes more braking in the ejecta and sends 
the reverse shock back through it more quickly. In con- 
trast, we note that Mfe is essentially constant when mass 
loading occurs by hydrodynamic ablation (Paper I). Since 
the mass loading is essentially saturated in both formu- 
lations at this stage, these differences arise from the fact 
that in this paper we also mass load in the shocked ejecta 



on n . For slightly higher mass loading rates (/' > 3.16) 
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Fig. 2. The ratio of evaporated or ablated mass to swept-up mass as a function of time (normalized to the time of shell 
formation cf. Eq. [SJ for mass loading by hydrodynamic ablation (top panels) or by conductive evaporation (bottom 
panels). t s { = 4 x 10 5 yr (n = 0.01) and 2 x 10 3 yr (n = 100). The diamonds mark the end of the QST stage, and 
almost overlap in b). 



(unlike in Paper I) , which has the effect of increasing its 
pressure. 

In Fig. |31 we show the radius at the end of the free- 
expansion stage, i?FE; as a function of n and /'. We 
derive an appropriate analytical approximation for the 
dependence of -Rfe on no and /' by following the same 
procedure as in Sec. 4.1 of Paper I. The agreement be- 
tween the simulations and the analytical approximation is 
generally good, although the latter systematically under- 
estimates i?FE at high values of /' (as was also found in 
Paper I, where an explanation for this is given). For future 
purposes, however, we are more concerned with accurate 
estimates for Rqst, the radius at the end of the Quasi- 
Sedov- Taylor phase. We note that the much higher values 



of / which we consider in this paper lead to a greater range 
in i?FE at a given uq. 

4.2. The Quasi-Sedov-Taylor (QST) Phase 

In Fig. and Table [21 we show the variation of -Rqst with 
no and /'. We find that the expression 

RQSTif, n ) « 0.75i? QST (/' = 0)/'-°- 05 10- fo ^Mno/a) ^ (12) 

where a = 6 x lQ-0-9(i-e- x )x ; b _ 017 (i ogl0 /' + 2.5) 11 , 
and x = logic/' + 0.5 achieves a good fit to the numerical 
results. We plot values of -Rqst using Eq.^Jin Fig. 

The ratio of the radius at the end of the Quasi-Sedov- 
Taylor phase to the radius at the same stage for the case 
where f — (^qst/-Rqst(/' = 0)) is shown in Fig.HUlfor 



Pittard et al. : Evolution of Mass Loaded SNR II. 



9 




Fig. 3. Mass fractions (left panels) and energy fractions (right panels) as functions of the remnant age for /' = 
0, 3.16, 31.6 (/ = 0, 10 4 , 10 5 ) and n = 0.01 {t st = 4 x 10 5 yr, E = 10 51 erg, M = 10 M©). In each of the panels on the 
left the solid line shows the mass fraction of evaporated material, the dashed line shows the ejecta mass fraction, and 
the dotted line shows the swept-up mass fraction. In each of the panels on the right the dashed line shows the kinetic 
energy fraction, and the dotted line shows the thermal energy fraction, both in terms of the current remnant energy. 
The solid line shows the total energy as a fraction of the initial remnant energy. Note that we subtract the thermal 
energy swept-up by the remnant - if this is included both the total energy and the thermal energy fraction increase at 
late times. The vertical lines mark the ends of the free expansion (FE) and Quasi-Sedov- Taylor (QST) phases. 
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Fig. 4. Remnant evolution of density, pressure, temperature, and velocity as a function of age for n = 0.01 and /' = 
(left panels) or /' = 3.16 (right panels). The profiles are at time t=3.54 x 10 4 yr (solid); t=8.59 x 10 4 yr (dotted); 
t=1.70 x 10 5 yr (dashed); t=3.39 x 10 5 yr (dot-dashed); t=5.41 x 10 5 yr (dot-dot-dot-dashed). (t s[ = 4 x 10 5 yr). 
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Fig. 5. Remnant averaged mass injection rate a), and integrated mass b), as a function of time for no = 0.01 cm~ 3 
and /' = 3.16 (solid), /' = 31.6 (dotted), and /' = 316 (dashed). Here we define the remnant volume as a sphere with 
radius equal to that of the blast wave. The slight increase in p at early time results from the fact that the post-shock 
region where mass loading occurs is above the saturation temperature, while the volume filling factor simultaneously 
increases as the reverse shock propagates back to the centre of the remnant. Note also that the maximum amount of 
mass that can be evaporated when no = 0.01 cm -3 is < 3000 M Q . 



all /' and no. Although there is variation with no for no < 
10, this ratio appears to become fairly insensitive to no for 
larger values of/'. We further note that the behaviour seen 
in Figs. 171 and 1 101 is qualitatively similar (which is also the 
case for ablative mass loading cf. Paper I). 

For an ambient density, n = 0.01, the approximation 
noted in Eq. 11 of Paper I, with substitution of / with /', 
i.e. 

i? Q ST « 0.9i? QST (/' = 0)(10v/7r - 09kW ', (13) 

is a good fit, while for high ambient densities (no ~ 100.0) 
an excellent approximation is 

logioi? QS T « -0.31og 10 /' - 0.5. (14) 

In the equations in this section it is again assumed that 
the remnants have an initial energy, E — 10 51 erg, and 
ejecta mass, M = 10 M Q . 

5. Conclusions 

We have extended work presented in Paper I on the evo- 
lution of mass loaded supernova remnants by considering 
mass loading by conductively driven evaporation of em- 
bedded clouds. Paper I confirmed that mass injection can 
strongly influence remnant evolution, and we also estab- 
lish here that the nature of the mass injection process 
is also important in this regard. This is due to the fact 
that mass loading through conductive evaporation is ex- 
tinguished at relatively early times. Hence, conductively 



driven mass loading does not appreciably alter the later 
stages of remnant evolution, when the remnant is domi- 
nated by swept up ambient gas. Therefore remnants that 
ablatively mass load are dominated by loaded mass and 
thermal energy at late times (Paper I), while those that 
conductively mass load are dominated by swept-up mass 
and kinetic energy. The greater dominance of loaded mass 
in the ablative case means that such remnants evolve more 
quickly, and reach all dynamical stages earlier. At a given 
age they tend to be both more massive and smaller than 
equivalent remnants which are conductively mass loaded. 

We are able to confirm some of the properties of 
conductively mass loaded remnants predicted from self- 
similar solutions, and in particular find that such rem- 
nants may display a thick-shell morphology (cf. the hy- 
drodynamic results presented in Cowie et al. H981l and the 
similarity solutions presented in Dyson & Hartquist .1987). 

In this work we have been particularly interested in the 
range of conductively mass loaded supernova remnants at 
the time at which they have radiated away half of their 
initial energy (see Table 01. It was noted in Sec. that /' 
may be dependent on no- This behaviour would pick-out 
a roughly diagonal line in Tabled and implies that the 
radius of remnants at the end of the quasi-Sedov- Taylor 
stage has less variance with no than would otherwise be 
the case. However, since /' is also dependent on the num- 
ber density of cold clouds we do not expect a particularly 
tight relationship. 
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Fig. 6. The fractional volume of the remnant above a temperature of 10 5 K (solid), 10 6 K (dotted), and 10 7 K (dashed) 
for no = 0.01 cm~ 3 and /' = 3.16 a) or /' = 31.6 b). At early times when the ejecta is being thermalized all of the 
shocked gas is hotter than 10 K and the curves are coincident in each case, and rising towards a value of 1.0. When 
the mass loading is stronger the remnant cools quicker. 



Simple approximations that fit the evolution of the 
range of supernova remnants which conductively mass 
load, and which are complementary to similar approxi- 
mations in Paper I, have also been found. In both works 
it is assumed that the remnant has an initial energy, 
E = 10 51 erg, and an ejecta mass, M = 10 Mq. We expect 
that the evolution of remnants which ablatively mass load 
will be fairly insensitive to the progenitor mass, because 
the majority of the mass loading occurs after the FE stage 
ends, at which point the swept up mass is about 6 times 
greater than the progenitor mass. In the conduction case 
it is not so clear what will happen, because most of the 
mass loading occurs in the early stages when the remnant 
is becoming thermalized. Furthermore, the mass loading 
may well be more dependent on our model assumptions 
(such as the saturation temperature, T sa t), than on the 
progenitor mass, at least in some regions of parameter 
space. With regards to the possibility of different explo- 
sion energies, we note that our solutions should scale in a 
similar way to the time of thin shell formation, t st cx 

£3/14 

(cf. Eq.©. 

Our range approximations will form the basis of fu- 
ture work to investigate galactic superwinds formed by 
the combination of many overlapping supernovae. 

Acknowledgements. We would like to thank Dr. R. Bandiera 
for constructive comments which led to clarification of our as- 
sumptions and generally improved the paper. JMP would also 
like to thank PPARC for the funding of a PDRA position. 
This work has made use of NASA's Astrophysics Data System 
Abstract Service. 



^0 00 F 1 — ■ I I 1 1 1 1 1 1 — I I I I 1 1 1| 1 — I I I ■ 1 1 ■ I 1 — I I I MIL 



1.00 r 



0.10 r 




"V. s 



0.01 1 i I i I i I i 

0.1 1.0 10.0 100.0 1000.0 

f 



Fig. 7. Ratio of the half energy time to the thin shell for- 
mation time (ti/2/tsf) as a function of the mass loading 
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Fig. 10. Ratio of the radius at the end of the Quasi-Sedov- 
Taylor phase to the radius at the same stage for the case 
where /' = (#qst/-Rqst(/' = 0)) against /' for all n : 
n = 0.01 (solid); n = 0.1 (dotted); n = 1.0 (dashed); 
n = 10.0 (dot-dashed); n — 100.0 (triple dot-dashed). 
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